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Summary  of  Objectives  and  Approach: 

Resonant  Ultrasound  Spectroscopy  (RUS)  has  been  used  successfully  to  determine  the  elastic 
properties  of  single  crystals  and  homogeneous  materials.  Our  objective  was  to  develop  RUS  to 
measure  the  average  elastic  properties  of  inhomogeneous  materials,  such  as  composites, 
concrete,  and  rocks.  Our  goal  was  to  proscribe  the  conditions  necessary  for  the  successful  use  of 
RUS  on  inhomogeneous  materials  by  studying  the  effects  of  size,  shape,  aspect  ratio,  and  sample 
Q  on  the  accuracy  of  RUS  measurements.  We  focused  on  rock  samples,  believing  that  the 
strategy  for  establishing  the  applicability  of  RUS  to  rocks  should  generahze  to  other 
inhomogeneous  materials.  Rocks  are  consolidated  materials,  not  easily  machined,  often 
anisotropic,  and  having  very  low  quality  factor,  Q  (or  high  attenuation). 

We  have  attempted  to  answer  the  following  question.  Under  what  conditions  is  RUS  a  useful 
tool  for  determining  die  moduli  of  macroscopic  (a  few  centimeters  on  a  side),  inhomogeneous 
samples?  We  concentrated  on  identifying  a  sample  geometry  that  will  maximize  success  with 
RUS.  The  work  consisted  of  numerical  modeling  of  sample  resonances  under  varying 
conditions,  and  empirical  testing  of  rock  samples.  Numerical  modeling  was  based  on  the 
Visscher  et  al.  [1991]  variational  technique,  and  standard  perturbation  techniques.  Some 
software  was  developed  at  UNR,  and  some  was  provided  by  Dynamic  Resonance  Systems 
(DRS),  a  commercial  provider  of  RUS  measurement  systems.  Empirical  testing  was  performed 
using  a  DRS  system.  The  experimental  sample  shapes  are  rectangular  parallelepipeds. 
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Summary  of  Major  Findings: 

1.  The  RUS  analysis  technique  assumes  sample  homogeneity.  By  a  back-of-the-envelope 
argument,  we  claim  that  inhomogeneous  samples  will  look  homogeneous  to  RUS  if  the 
scale  of  the  inhomogeneity  is  no  more  than  1/5  the  smallest  sample  dimension.  Thus 
average  moduli  can  be  accurately  determined  when  the  scale  of  sample  inhomogeneity  is 
small  compared  to  the  sample  scale.  Empirical  testing  on  inhomogeneous  granite 
samples  of  various  sizes  confirmed  this  concept. 

2.  The  RUS  analysis  technique  assumes  a  perfect  parallelepiped  sample.  Using  perturbation 
theory,  we  argue  that  for  inhomogeneous  samples,  a  nick  in  a  comer  of  the  sample,  or 
two  sample  sides  that  are  not  perfectly  parallel,  will  contribute  negligibly  to  error  in  the 
measured  moduli.  Empirical  results  from  crudely  cut  samples  and  samples  cut  precisely 
are  essentially  indistinguishable  for  the  same  rock  type.  Small  sample  shape 
imperfections  did  not  cause  large  errors  in  fits  for  moduli. 

3.  Inhomogeneous  samples  have  high  attenuation,  or  low  quality  factor  Q.  Low  Q  causes 
resonance  peaks  to  be  short  and  broad.  This  means  that  resonance  peaks  may  overlap, 
and  may  be  difficult  to  pick  out  from  background  noise  (especially  as  the  frequency 
increases).  To  maximize  the  information  gained  from  a  minimum  number  of  peaks,  we 
modified  sample  aspect  ratio.  We  modeled  the  spread  in  the  first  14  resonance 
frequencies  as  a  function  of  aspect  ratio,  and  die  contribution  of  the  two  isotropic  moduli 
to  the  first  resonances  as  a  function  of  aspect  ratio.  In  both  cases,  we  found  fliat  optimum 
results  can  be  obtained  with  sample  aspect  ratios  of  4/1  (largest  to  smallest  side). 
Empirical  testing  of  black  gabbro  samples  with  varying  aspect  ratios,  show  that  the  error 
in  the  fits  for  moduli  decreases  as  the  aspect  ratio  is  increased. 

4.  In  theory,  the  RUS  technique  should  be  capable  of  identifying  anisotropy  in  an 
experimental  sample.  In  practice,  inhomogeneity  and  compensations  made  to  deal  with 
inhomogeneity  may  mask  clues  that  a  sample  is  anisotropic.  In  other  words,  if  we 
assume  isotropy  for  an  anisotropic  sample,  the  RUS  fitting  errors  should  be  large  and 
obvious  for  a  perfect  sample,  but  may  not  be  for  inhomogeneous  samples.  We  modeled 
the  error  induced  in  an  isotropic  fit,  as  a  function  of  the  strength  of  hexagonal  anisotropy. 
We  conclude  that  anisotropy  will  be  very  difficult  to  detect  in  high  aspect  ratio  samples. 
In  contrast,  anisotropy  in  symmetric  samples  will  cause  degenerate  resonances  to  split, 
and  should  be  clearly  distinguishable.  However,  low  Q  may  cause  these  split  peaks  to 
overlap,  and  thus  make  the  distinction  impossible. 

Summary  of  Conclusions: 

Numerical  modeling  and  empirical  testing  indicate  that  RUS  is  a  viable  technique  for. 
characterizing  the  average  isotropic  elastic  moduli  of  inhomogeneous  materials,  although  larger 
RMS  errors  can  be  expected  than  for  single  crystal  materials.  Success  with  RUS  can  be 
optimized  by  ensuring  that  the  sample  size  is  large  compared  to  the  scale  of  inhomogeneity,  and 
by  using  a  high  aspect  ratio  parallelepiped  sample  (ratio  of  largest  side  to  smallest  side 
approximately  4).  Hexagonal  anisotropy  and  lower  symmetries  may  be  difficult  to  detect,  and 
even  harder  to  quantify,  except  for  very  low  attenuation  materials. 


Detailed  Report: 

See  attached  Master’s  Thesis,  Determination  of  Elastic  Moduli  of  Rock  Samples  Using  Resonant 
Ultrasound  Spectroscopy,  by  TJ  Ulrich,  Department  of  Physics,  University  of  Nevada,  Reno, 
December  2000. 

Navy  Relevance: 

The  development  of  RUS  as  a  tool  to  understand  thermoacoustic  properties  of  crystalline 
materials  is  supported  primarily  by  the  Physical  Acoustics  Program  at  ONR.  Our  goal  is  to 
develop  RUS  as  a  tool  to  understand  the  elastic  and  acoustic  properties  of  materials  that 
constitute  physical  infrastructure,  such  as  concrete,  rock,  laminates,  and  polymer  composites. 
This  corresponds  to  ONR’s  Physical  Acoustics  Program  objective  of  measuring  and  correlating 
physical  properties  of  materials  using  acoustic  methods. 
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Continuing  Program: 

During  the  past  year,  we  performed  a  feasibility  study  showing  that  with  the  proper  choice  of 
sample  geometry,  RUS  can  be  used  to  measure  elastic  moduli  of  inhomogeneous  materials 
quickly  and  accurately.  The  RUS  measurement  system  can  be  used  to  monitor  the  elastic  state 
of  an  inhomogeneous  sample.  In  the  next  year,  we  have  two  goals:  1)  to  broaden  the  range  of 
inhomogeneous  samples  to  which  RUS  can  be  applied  (we  are  particularly  interested  in 
cylindrical  samples),  and  2)  to  develop  RUS  as  a  tool  to  probe  the  response  of  rocks  and  concrete 
to  quasistatic  and  dynamic  perturbations  (we  are  particularly  interested  in  dependence  on 
temperature  and  vapor  pressure).  An  experimental  system  will  be  developed  to  allow  control  of 
temperature  (first)  and  vapor  pressure  (second),  and  a  series  of  experiments  that  follow  the 
evolution  of  the  elastic  state  in  response  to  these  variables  will  be  performed.  We  believe  the 
outcome  of  these  experiments  will  provide  important  new  insights  regarding  the  microscopic 
mechanism  causing  the  nonlinear  elastic  behavior  displayed  by  rocks  and  concrete. 
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Abstract 

Resonant  Ultrasound  Spectroscopy  (RUS)  is  a  method  whereby  the  elastic 
tensor  of  a  sample  is  extracted  from  measured  resonance  frequencies.  RUS  has  been 
used  successfully  to  determine  the  elastic  properties  of  single  crystals  and  homoge¬ 
neous  samples.  In  this  thesis,  I  study  the  applicability  of  RUS  to  macroscopic  samples 
of  inhomogeneous  materials,  specifically  rock.  I  have  paid  particular  attention  to  sev¬ 
eral  issues:  appropriate  scale  of  inhomogeneity,  sample  shape  imprecision,  effects  of 
low  Q,  extraction  of  siiflScient  information  on  multiple  moduli,  and  recognition  of 
anisotropy.  Using  modeling  and  empirical  testing,  I  have  established  bounds  on  the 
applicability  of  RUS  to  inhomogeneous  materials. 
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“Six  months  in  the  lab  can  save  you  a  day  in  the  library.” 


-  old  adage 
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Chapter  1 

Introduction 


Resonant  ultrasound  spectroscopy  (RUS)  is  a  highly  sensitive  technique  for 
determining  the  elastic  tensor  and  related  properties  of  single  crystal  samples  [1].  The 
first  ten  or  more  resonance  frequencies  of  a  sample  with  free  boundaries  are  measured. 
These  resonance  frequencies  are  the  input  to  an  iterative  inversion  algorithm  that 
finds  the  best  match  between  the  data  and  a  set  of  resonance  frequencies  generated 
from  a  model.  The  variational  parameters  in  the  inversion  are  the  components  of 
the  elastic  tensor  describing  the  sample.  Thus  RUS  is  a  method  whereby  the  elastic 
tensor  of  a  sample  is  extracted  from  measured  resonance  frequencies.  For  example, 
RUS  has  been  used  to  determine  the  elastic  tensor  and  crystallographic  orientation  of 
tantalum  [2],  to  measure  the  elastic  properties  and  infer  thermodynamic  properties  of 
CaO  at  temperatures  representative  of  the  earth’s  interior  [3],  and  to  study  the  elastic 
constants,  microstructure,  and  superconducting  phase  transition  in  La2_xSr3;Cu04  [4]. 
The  elastic  tensor  of  single  crystals  and  minerals  reflects  the  angstrom  scale  structure 
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of  these  systems.  From  knowledge  of  the  elastic  tensor,  the  symmetry  of  a  crystal  and 
features  of  the  phonon  and  electron  spectra  can  be  deduced.  When  RUS  is  conducted 
in  a  temperature  reservoir  that  can  be  controlled,  it  is  possible  to  develop  a  complete 
thermodynamic  description  of  the  sample. 

The  success  of  RUS  derives  from  the  sensitivity  of  the  resonant  modes  of  a  sys¬ 
tem  to  the  elastic  structure  of  the  system  [5].  The  elastic  structure  affecting  resonance 
frequencies  has  three  components;  the  geometry  of  the  sample;  the  homogeneity  of 
the  sample;  and  the  elastic  tensor  of  the  sample,  including  symmetry  and  orientation. 
Given  a  perfectly  homogeneous  sample  with  a  precise  geometry,  the  elastic  tensor  can 
be  derived  to  a  very  high  degree  of  accuracy. 

In  this  paper,  we  explore  the  limits  of  applicability  of  RUS  to  macroscopic,  in¬ 
homogeneous  samples  of  rock.  Rocks  are  consolidated  materials,  typically  assembled 
from  aggregates  of  mesoscopic  sized  pieces  of  atomically  uniform  materials.  They  are 
not  easily  machined  to  precise  shapes.  Their  microscopic  symmetry  is  homogenized 
by  the  process  of  their  macroscopic  assembly.  However,  they  may  have  macroscopic 
symmetry  of  great  importance.  My  goal  is  to  proscribe  the  conditions  necessary  for 
the  successful  use  of  RUS  on  inhomogeneous  materials.  These  conditions  include  con¬ 
straints  on  sample  preparation  and  contraints  on  the  body  of  reasonable  questions 
that  can  be  answered  with  RUS. 

Historically,  the  elastic  moduli  of  macroscopic  inhomogeneous  materials  have 
been  measured  using  mechanical  testing  [6]  or  ultrasonic  time-of-flight  techniques  [7]. 
Mechanical  testing,  wherein  the  stress-strain  relationship  is  measured  between  am- 
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bient  conditions  and  failure,  is  a  widely  accepted  method  for  measuring  material 
strength  and  toughness.  A  modulus  is  determined  from  the  stress-strain  slope  at 
low  strain.  The  greatest  disadvantage  to  using  mechanical  testing  to  determine  lin¬ 
ear  elastic  moduli  is  that  the  sample  is  destroyed  or  altered  as  a  result  of  the  test. 
Thus  results  cannot  be  confirmed  for  a  given  sample,  and  only  one  elastic  constant 
can  be  determined  for  each  sample.  Another  disadvantage  is  that  mechanical  tests 
are  inherently  high  amplitude  tests.  Our  primary  interest  is  in  dynamic  moduli  that 
are  directly  related  to  acoustic  velocities.  In  general,  dynamic  moduli  measured  by 
time-of-flight  techniques  are  higher  than  moduli  infered  from  mechanical  testing. 

Time-of-flight  determinations  of  elastic  constants  are  actually  acoustic  velocity 
measurements.  The  time  delay  t  of  an  ultrasonic  pulse  across  the  sample  is  measured. 
Given  the  length  I  and  density  p  of  the  sample,  the  wave  velocities  and  elastic  con¬ 
stants  can  be  determined  from  v  =  //<,  cn  =  pv^,  and  C44  =  where  Vc  and  Vg 
are  the  compressional  and  shear  wave  velocities  (Table  1.1  provides  a  relationship 
between  cn,  C44  and  the  standard  engineering  elastic  constants).  This  technique  is 
path  dependent  and  therefore  requires  multiple  measurements  in  different  locations 
and  directions  to  determine  an  average  velocity  if  the  material  is  inhomogeneous.  To 
determine  both  compressional  and  shear  velocities,  transducers  that  produce  purely 
compressional  and  purely  shear  plane  waves  are  desirable.  Finally,  the  transducers 
must  be  bonded  to  the  sample  reproducibly  for  all  measurements. 

Resonant  ultrasound  spectroscopy  is  a  potential  alternative  technique  for  de¬ 
termining  elastic  moduli  of  inhomogeneous  materials.  Typically,  a  rectangular  par- 
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allelepiped  sample  is  placed  between  two  piezoelectric  transducers,  a  source  and  a 
detector.  The  sample  is  driven  at  constant  voltage  as  the  frequency  is  swept  through 
multiple  resonances.  The  measured  resonance  frequencies  are  the  input  to  an  iter¬ 
ative  inversion  algorithm  that  finds  the  best  match  between  the  data  and  a  set  of 
resonances  generated  from  a  model.  The  experimental  component  of  RUS  is  straight¬ 
forward.  The  analysis  component  is  sophisticated,  but  can  be  performed  more  and 
more  rapidly  as  computational  power  increases.  Chapter  2  will  examine  the  numerical 
process  involved  in  the  RUS  analysis.  In  Chap.  3,  several  issues  pertaining  to  using 
RUS  on  inhomogeneous  samples  are  discussed.  In  most  cases,  numerical  modeling 
was  used  to  explore  ways  to  optimize  experimental  chances  for  success.  In  Chap.  4, 
the  results  of  RUS  experiments  on  a  variety  of  samples  are  displayed  and  discussed 
as  well  as  a  summary  of  our  findings,  describing  the  bounds  on  RUS  applicability  to 
inhomogeneous  materials  found  empirically  and  through  modeling. 
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Chapter  2 


Variational  Calculation  of 
Elastic  Resonances  and 
Levenberg-Marquardt  Inversion 


The  Resonant  Ultrasound  Spectroscopy  (RUS)  technique  has  two  components: 
the  experiment  to  collect  resonance  frequencies,  and  the  analysis  to  extract  the  elastic 
moduli.  The  bulk  of  this  thesis  is  a  description  of  the  experimental  application  of 
RUS,  and  the  numerical  modeling  used  to  optimize  the  experimental  results.  In 
this  chapter  the  mathematical  and  theoretical  basis  of  the  analysis  component  will  be 
described  for  two  physical  systems:  a  one  dimensional  string,  and  a  three  dimensional 
elastic  solid.  The  example  in  one  dimension  is  used  to  illustrate  the  principles  of  the 
three  dimensional  calculation. 

While  elastic  resonances  in  one  dimension  are  relatively  easy  to  calculate  ana- 


( 

lytically,  the  resonances  of  an  arbitrary  three  dimensional  object  with  arbitrary  sym¬ 
metry  are  extremely  difficult  to  calculate.  The  analysis  built  into  the  RUS  technique 
uses  a  variational  scheme  to  calculate  approximate  resonances  [8].  This  approximation 
technique  decreases  the  calculation  time  considerably  and  can  be  applied  to  arbitrary 
shape  and  symmetry  in  a  straightforward  manner.  In  addition,  extremely  accurate 
results  can  be  obtained  when  an  appropriate  number  of  basis  functions  is  employed. 


2.1  Variational  Calculation  of  Elastic  Resonances 
in  One  Dimension 


Hamilton’s  variational  principle  is  used  to  calculate  the  resonance  frequencies  of  an 
arbitrary  body.  We  illustrate  this  calculation  in  one  dimension.  The  Lagrangian  of  a 
string  of  length  I  is 

L=  j{T-  U)  dx  ,  (2.1) 

i 

where  T  is  the  kinetic  energy  and  U  is  the  potential  energy  per  unit  length  of  the 
string.  In  frequency  space 

T  =  (2.2) 


and 


(2.3) 
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where  p  is  the  mass  density,  r  is  the  string  tension,  u  is  the  frequency  of  motion,  and 
u  is  the  transverse  displacement  at  x.  Thus  the  Lagrangian  is 

L  =  I  -  ir  I  dx  .  (2.4) 

I  I  ^  J 

Varying  our  original  Lagrangian  by  6L,  Eq.  (2.4)  becomes 

L  +  6L=I  {u  +  6uf  -  ]  cte-  (2-5) 

To  first  order  in  the  variational  parameter  du, 


Integrating  by  parts  we  find 

i 

(2.7) 

0 

To  minimize  the  Lagrangian  we  set  6L  =  0.  For  arbitrary  bu,  Hamilton’s 
principle  is  satisfied  when  two  conditions  are  met: 

First  Condition: 

I  » 

=  0  (2.8) 

0 

Second  Condition: 

puj%.  +  r  =  0  (2.9) 

The  first  condition  is  equivalent  to  the  free  boundary  condition  at  each  end  of 
the  string.  The  second  is  the  one  dimensional  wave  equation.  Thus  the  displacement 
function  u{x)  that  minimizes  the  Lagrangian,  solves  the  elastic  wave  equation  for  a 
string  with  free  ends,  as  illustrated  in  Fig.  2.1. 
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Applying  the  results  of  the  previous  paragraph,  we  choose  a  convenient  form  for 
v{x)  that  includes  variational  parameters  that  will  allow  us  to  minimize  L.  Following 
Visscher  et  al.  [9]  we  choose  u{x)  to  be  a  polynomial, 


N 


n=0 


Substituting  u{x)  in  Eq.  (2.4)  we  have 

N  N 


r  1  ^  ^  1  ^  dx”  ^ 

L  =  [  dx  -pj^  ^  a„.T”  ^  arr,.x^  -  -T  an 

I  n=0  rn=0  n=0  m=0 


dx^ 

dx 


Simplifying, 


1  AT  1  iv 

L  —  ^  ^  ajjajjiEnm  q  ^  ]  andf^n  m  rn,ni.  > 

^  n,m=0 


N 


n, 777=0 


where 


En.m.  — 


pln+m.+l 

n  +  m  +  1' 


and 


r„„  =  < 


nm 


in+n<  1  for  n  &  m  ^  0 


(2.10) 


(2.11) 


(2.12) 


(2.13) 


(2.14) 


r?+m- 

0  for  n  or  m  =  0 

To  minimize  the  Lagrangian,  we  set  dL/daj  =  0.  This  leads  to  N  linear 

equations  that  can  be  expressed  in  matrix  form. 


(a;2E-r)a  =  0. 


(2.15) 


We  now  have  an  eigenvector /eigenvalue  problem.  Setting  the  determinant  |w^E— Fj  — 
0,  we  find  A^+ 1  values  of  a;^.  For  each  frequency  eigenvalue  there  is  an  eigenvector 
a  that  contains  the  values  of  the  coefficients  of  the  displacement  polynomial  u{x). 

As  a  concrete  example,  take  N  =  2.  Then  u{x)  =  oq  +  aix  +  a2X^.  Let 
/  =  p  =  r  =  1  for  simplicity.  Then  the  eigenvalue/eigenvector  equation  reduces  to 
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to 

3 

1  ,  ,2 

ia;2-l 

ia>2-l 

1  ,  ,2  4 

5  3 

-  r  = 


Setting  |u;^E  -  r|  =  0,  we  find  the  eigenvalue/eigenvector  pairs: 

/1\ 

0 

\0/ 

(  1  \ 


u}q  —  0 ,  ao  — 


U! 


^  =  12 ,  ai  = 


ul  =  60 ,  a2  = 


-2 

V  0  / 

/  1  \ 


V  6  / 


The  nontrivial  eigenvalue/eigenvector  pairs  give  us  the  two  solutions 

ui{x)  =  1  —  2x, 


(2.16) 


(2.17) 


(2.18) 


(2.19) 


(2.20) 


and 

U2{^')  =  1  —  6x  +  6x‘^  .  (2.21) 

Comparing  these  numerical  solutions  to  the  anal3d;ical  solutions 

=  (^)  =  cos{ujx),  (2.22) 

(Fig.  2.2),  we  can  see  that  the  shape  is  not  identical.  As  N  increases,  the  nature  of 
the  variational  technique  is  such  that  the  frequencies  from  the  variational  solution 
converge  rapidly  to  those  of  the  exact  solution,  despite  the  inaccuracies  in  the  modes 
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shapes  [10].  Notice  that  the  first  frequency,  wg,  is  zero.  This  corresponds  to  the  trans¬ 
lation  of  the  string  with  no  vibration.  As  an  experiment  only  measures  vibrational 
modes,  we  can  neglect  the  translational  mode.  It  is  apparent  now  that  assuming  a 
polynomial  solution  of  order  A  =  2  is  very  restrictive,  limiting  us  to  N  frequencies 
that  can  be  calculated  and  compared  to  the  experiment.  Note  that  an  N—  order 
polynomial  solution  for  u{x)  has  a  maximum  of  N  nodes. 


2.2  Three  Dimensional  Model 


Figure  2.3  shows  some  exagerated  mode  shapes  for  a  parallelepiped.  Comparing 
these  to  the  modes  shapes  shown  in  Fig.  2.2,  it  is  immediately  obvious  that  the 
three  dimensional  model  will  be  much  more  complicated.  Hamilton’s  variational 
principle  can  be  used  to  find  the  resonance  frequencies  of  a  3-D  elastic  solid  in  much 
the  same  way  as  for  the  1-D  case.  However,  the  complexity  of  the  problem  rises 
dramatically,  especially  for  non-isotropic  systems.  Here  we  state  the  results  and 
equations  equivalent  to  those  described  in  Sec.  2.1. 

We  choose  a  three  dimensional  polynomial  for  the  variational  solution. 


Ui  =  aa  ,  (2.23) 

A 

where 

=  (2.24) 

N  >  I  +  m.  +  n  ,  (2.25) 


and  A  is  an  index  to  distiguish  between  sets  of  l,m,n.  The  Lagrangian  of  Eq.  (2.1) 


Figure  2.3;  Examples  of  three  dimensional  mode  shapes. 


for  a  three  dimensional  system  then  becomes 


V-  c  2  a.  a.  1  V- 

^  buifKjj^aixai'x^x^x'  -  x  Cijkiaixai/y 


dxj  dxi 


dV,  (2.26) 


where  Cjjki  is  defined  from  the  generalized  Hooke’s  law  [11] 


<Xij  —  Cijki  eki, 


(2.27) 


and 


—  ^jikl  ('k.lij  ^ijlkj 


(2.28) 


by  symmetry.  Thus,  only  21  of  the  81  components  are  independent,  and  Cijki  is 
commonly  rewritten  as  a  36-component  second  order  tensor  Cij.  Also  by  symmetry, 
ffij  and  Cki  have  only  six  independent  components,  and  are  commonly  rewritten  as 


6-component  first  order  tensors. 
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Defining 

Exix'i'  =  6ii'  j  ^xP^x'  dV  (2.29) 

Jv 


and 


^  V-  /"  JT/ 

^MX'i'  -^Cijer  1^^  g  Q,  dV, 
J.J  ■'  •' 


we  can  write  the  lagrangian  in  tensor  form, 


L  =  -  ^a^Fa, 


1 


(2.30) 


(2.31) 


analogous  to  Eq.  (2.12)  for  the  1-D  case,  but  more  difficult  to  perform  as  a  is  now  a 
rank  3  tensor  with  elements.  This  then  leads  to  an  eigenvalue  equation 


a;2Ea  =  Ta.  (2.32) 

Solving  this  system  as  >  oo  would  give  an  exact  solution,  but  is  analytically 
impossible.  Thus  a  numerical  solution  must  be  found.  To  allow  for  a  realistic  com¬ 
putational  time  we  have  limited  our  numerical  solution  to  a  polynomial  solution  of 
the  order  A^  =  10. 

For  any  nontrivial  N,  the  diagonalization  of  Eq.  (2.32)  is  performed  numer¬ 
ically.  As  in  the  one  dimensional  case,  the  variationally  determined  frequencies  in¬ 
crease  in  accuracy  as  the  parameter  N  increases.  However,  the  computational  time 
and  memory  requirements  also  increase  with  N.  We  have  limited  the  order  of  our 
three-dimensional  polynomial  to  N  —  10.  Figure  2.4  [12]  shows  a  comparison  of 
the  number  of  resonance  frequencies  that  may  be  determined  accurately  for  A^  =  10 
and  N  =  12.  The  sample  is  SrTiOs  and  of  cubic  symmetry  (3  independent  elastic 
constants).  The  numerically  derived  frequencies  are  compared  to  the  experimentally 


dN/df  (lines/0.5MHz) 
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SrTI03  Sample  Cl ,  r-  3kHz 
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Figure  2.4:  Comparison  of  numerical  method  for  number  of  frequencies  predicted  by  N  = 
10  and  AT  =  12  to  the  actual  number  of  observed  frequencies.  [12] 
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derived  frequencies  and  to  the  standard  p  result.  Note  that  up  to  50  frequencies  may 
be  fit  with  the  choice  N  =  10.  Since,  in  general,  we  cannot  distinguish  more  than  20 
resonances  from  our  rock  samples,  we  are  comfortable  with  N  =  10. 

2.3  Levenberg-Marquardt  Minimization 

In  Sec.  2.2  we  describe  a  forward  calculation  of  resonance  frequencies.  After  the 
resonance  frequencies  have  been  calculated  for  a  given  initial  guess  of  Cijki  (r  from 
Sec.  2.1),  the  RUS  inversion  algorithm  compares  the  calculated  values  to  those  mea¬ 
sured  experimentally.  If  the  calculated  frequencies  do  not  match  the  experimental 
values,  a  new  elastic  tensor  is  chosen  using  a  Levenberg-Marquardt  minimization 
of  where  is  defined  as  the  rms  error  between  the  calculated  and  measured 
frequencies, 

M 

=  Ufi"  -  ffW  (2-33) 

7=1 

M  is  the  number  of  frequencies  used  for  the  fit,  is  the  i**’  experimentally  measured 
frequency,  and  ff  is  the  calculated  frequency.  These  three  steps:  calculation  of 
resonances  from  an  elastic  tensor,  comparison  of  calculated  and  experimental  reso¬ 
nances,  and  a  new  choice  of  elastic  tensor,  are  repeated  until  a  fit  is  achieved.  What 
is  defined  as  a  good  match  or  fit  will  be  discussed  in  Ch.  3.  Here  we  will  discuss  some 
of  the  details  of  the  minimization  technique. 

The  Levenberg-Marquardt  minimization  technique  [13]  combines  a  steepest 
decents  minimization  with  the  Newtonian  minimization  method.  This  allows  a  rapid 
convergence  upon  the  absolute  minimum  of  the  surface.  The  steepest  descents 
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method  is  used  far  from  the  minimum  in  The  Levenberg-Marquardt  algorithm 
switches  to  the  Newtonian  method  continuously  as  the  minimum  is  approached, 
is  a  dimensionless  quantity,  and  is  looked  upon  throughout  this  thesis  as  a  surface  in 
Cii,  C44  space.  This  analogy,  however,  does  not  hold  as  the  modeled  material  deviates 
from  isotropy,  i.e.,  as  more  independent  elastic  constants  are  introduced. 

The  general  procedure  for  using  the  Visscher  et  al  [9]  variational  technique 
and  the  Levenberg-Marquardt  method  to  determine  elastic  constants  of  a  material  is 
as  follows: 

•  Guess  initial  values  of  elastic  constants  {C guess)- 

•  Compute  resonance  frequencies  for  a  sample  characterized  by  (Cg„ess)  and  the 
geometry  of  the  experimental  sample. 

•  Compute  guess,  Cexp)^ox  the  model  frequencies. 

•  Increment  Cguess  and  repeat  the  procedure  above  until  x^iCguess,  Cexp)  is  min¬ 
imized. 

In  this  procedure,  the  increment  size  of  Cguess  can  vary,  and  depends  on  the  value  of 
Cexp)-  This  allows  for  a  rapid  and  accurate  convergence  upon  the  minimum 
even  when  starting  at  a  point  far  from  the  absolute  minimum.  It  should  be  noted, 
however,  that  due  to  the  potential  complexity  of  the  x^  surface,  relative  or  false 
minima  may  exist.  For  this  reason  extreme  differences  in  the  initial  guess  and  the 
absolute  minimum,  may  cause  a  false  minimum  to  be  found  and  result  in  erroneous 
results;  thus  the  more  educated  we  can  be  when  making  our  initial  guess  at  the  elastic 
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constants,  the  more  certain  we  can  be  that  we  have  found  the  absolute  minimum  [13]. 

Another  important  feature  of  the  Levenberg-Marquardt  method,  is  that  it 
provides  information  about  the  surrounding  surface  of  about  the  mininum.  That 
is  to  say,  that  the  curvature  or  sharpness  of  the  minimum  is  known  when  the  minimum 
is  found.  This  sharpness  can  differ  in  orders  of  magnitude  from  sample  to  sample, 
which  allows  us  to  not  only  find  the  elastic  constants,  but  to  assess  the  degree  of 
certainty  in  our  results.  A  very  shallow  minimum  implies  there  is  room  for  error, 
while  a  sharp  minimum  makes  us  much  more  confident  in  our  results.  This  is  an 
important  key  when  evaluating  the  quality  of  our  experimental  results,  and  will  be 


discussed  more  in  Sec.  3.6. 
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Chapter  3 


Modeling  and  Experimental 
Development 


In  this  section  we  will  apply  the  methods  described  in  Chap.  2  to  model  and 
analyze  experiments  performed  on  macroscopic  samples  of  inhomogeneous  rock.  The 
assumptions  inherent  in  the  analysis  described  in  Chap.  2  will  be  discussed,  as  well 
as  ways  in  which  to  maximize  the  success  of  RUS  on  samples  whose  properties  do  not 
satisfy  these  assumptions. 

The  experimental  system  is  shown  in  Fig.  3.1.  The  sample,  a  rectangular 
parallelepiped,  is  held  delicately  between  two  piezoelectric  transducers,  one  acting  as 
the  source  and  the  other  as  a  detector.  The  source  is  driven  at  constant  voltage  and 
the  frequency  is  swept  through  the  low-lying  resonances.  The  measured  resonance 
frequencies  are  the  input  to  the  iterative  inversion  algorithm  discussed  in  Chap.  2. 
The  algorithm  finds  the  best  match  between  the  data  (typically  10  to  20  resonance 
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Figure  3.1:  Diagram  of  experimental  set-up.  The  control  hardware  and  software  aie  com¬ 
mercial  packages  from  Dynamic  Resonance  Systems,  Inc. 


frequencies)  and  a  set  of  resonances  generated  from  a  model  in  which  the  parameters 

are  the  components  of  the  elastic  tensor. 

The  numerical  inversion  is  performed  by  software  provided  by  Dynamic  Res¬ 
onance  Systems  (DRS),  a  commercial  provider  of  RUS  measurement  systems,  and  is 
based  on  the  Visscher  et  al.  [9]  variational  technique.  The  analysis  software  employs  a 
model  of  the  resonant  elastic  system  that  (a)  has  free  boundaries,  (b)  is  spatially  ho¬ 
mogeneous,  and  (c)  is  a  rectangular  parallelepiped.  Departiire  of  a  sample  from  these 
.three  conditions  can  introduce  shifts  in  the  experimental  resonances,  and  thus  will 
introduce  error  in  the  derived  elastic  modtili.  Conditions  (a),  (b),  and  (c),  how  well 
our  samples  conform  to  them,  and  estimates  of  the  induced  error  due  to  nonconfor¬ 
mity  are  discussed  below.  In  Appendix  A  we  sketch  a  perturbation  theory  calculation 
that  we  use  to  assess  the  consequences  of  shape  departures  from  a  parallelpiped. 

Rock  and  similar  type  samples,  e.g.,  concrete,  have  relatively  high  acoustic 
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attenuation,  or  low  Q.  Thus  several  very  practical  issues  arise,  (d)  How  do  we  cause 
resonance  peaks  to  separate  enough  to  be  distinguishable?  (e)  How  do  we  acquire 
the  most  information  out  of  the  first  measurable  resonance  peaks,  i.e.,  how  do  we 
maximize  the  dependence  of  those  modes  on  the  elastic  constants  of  interest?  (f) 
To  what  extent  can  anisotropy  be  recognized  and  determined  for  an  inhomogeneous 
sample?  I  will  show  here,  that  by  altering  the  sample  geometry,  while  maintaining  a 
rectangular  parallelepiped  shape,  these  practical  issues  can  be  addressed. 


3.1  Free  boundaries 

The  variational  technique  used  in  the  analysis  [9]  is  based  on  the  recognition  that 
the  displacement  vectors  satisfying  the  elastic  wave  equation  with  free  boundaries  on 
the  sample  surface,  also  make  the  elastic  Lagrangian  of  the  sample  stationary.  To 
approximate  free  boundaries  in  the  experiment,  the  source  and  detector  are  most 
often  placed  at  vertices  of  the  parallelpiped,  delicately  supporting  the  sample.  The 
sample  is  nearly  free  standing.  Holding  the  sample  at  vertices  has  the  further  advan¬ 
tage  of  keeping  the  transducers  away  from  the  expected  node  lines  of  the  resonant 
modes.  When  the  full  resonance  spectrum  is  complicated,  transducers  may  be  placed 
purposely  at  expected  nodes,  such  as  the  center  of  a  face,  to  temporarily  simplify  the 
spectrum.  Using  transducers  for  support  limits  the  sample  size.  Our  transducers  are 
PZT-5  piezoelectric  pinducers.  We  have  limited  our  samples  to  less  than  50  cm^,  and 
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3.2  Inbomogeneity 
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3.3  Sample  geometry,  the  figure  of  the  sample 

Samples  of  consolidated  materials  are  difficult  to  machine  without  chipping,  and  often 
do  not  have  perfectly  parallel  sides.  How  ideal  must  the  figure  of  a  sample  be?  This 
question  can  be  examined  using  the  perturbation  treatment  of  the  elasticity  problem 
sketched  in  Appendix  A.  To  simulate  the  effect  of  an  error  in  the  figure  of  a  sample,  a 
localized  mass  is  carried  around  the  perimeter  of  the  sample,  and  the  frequency  shift 
caused  by  this  mass  is  calculated.  The  frequency  shift  for  mode  n  is  given  by 

^  «  2^  =  {Un\sp\u„),  (3.1) 

fn 

where  and  Un  are  the  eigenmodes  and  eigenfrequencies  of  a  perfectly  shaped  sam¬ 
ple,  and^p  is  the  localized  mass  perturbation  being  carried  around  the  sample.  In 
Fig.  3.2,  the  average  frequency  shift  of  the  first  20  resonances, 


where  N=20,  is  shown  as  a  function  of  the  perturbation  placement.  The  perturbation 
is  carried  along  the  sample  edge  and  into  the  sample  interior  as  shown  in  the  inset 
in  the  figure.  When  the  perturbation  is  at  an  interior  point  it  is  essentially  a  1% 
mass  distortion,  when  it  is  along  the  perimeter  it  is  a  1%  distortion  of  the  figure. 
Distortions  in  the  figure  of  the  sample  are  much  more  important  than  equivalent 
mass  distortions  in  the  sample  interior.  A  1%  chip  out  of  the  corner  of  a  sample  can 
produce  a  1%  change  in  the  frequency.  A  1%  mass  distortion  at  the  sample  center 
produces  less  than  0.2%  change  in  frequency. 

Given  the  number  of  other  contributors  to  error  in  RUS  measurements  on 
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Figure  3.2:  Frequency  shift  as  a  function  of  the  perturbation  placement. 

consolidatefl  materials,  these  contributions  to  error  are  rather  small.  This  conclusion 
was  confirmed  empirically  by  making  RXJS  measurements  on  samples  before  and  after 
chipping,  and  on  various  samples  of  tlie  same  size. 


3.4  Distinct  resonance  peaks 

Consolidated  materials  are  often  characterized  by  a  low  quality  factor  Q,  i.e.,  a  high 
attenuation.  At  fixed  amplitude  low  Q  materials  have  fewer  observable  resonance 
frequencies  than  high  Q  materials.  Additionally  the  broader  resonance  peaks  of  low  Q 
materials  overlap  nearby  peaks  and  complicate  peak  picking.  However,  the  geometry 
of  a  sample  sets  the  frequency  difference  between  peaks.  For  example,  a  sample  that 
is  a  cube  of  an  isotropic  material  has  a  three-fold  degeneracy  in  all  of  its  resonance 
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frequencies.  Thus,  we  can  use  geometry  to  minimize  peak  overlap  due  to  a  low  Q. 

In  Fig.  3.3,  calculated  resonance  freqtiencies  are  plotted  for  a  parallelepiped 
sample  as  a  function  of  the  aspect  ratio,  c/a.  The  volume  of  the  sample  is  fixed, 
a  X  6  X  c  =  4.8  cn?.^,  and  b  =  1.1a.  Aspect  ratios  greater  than  one  correspond  to 
rod-like  samples  and  are  characterized  by  the  number  c/a  in  Fig.  3.3.  Aspect  ratios 
less  than  one  correspond  to  plate-like  samples  and  are  characterized  by  the  number 
-a/c  in  Fig.  3.3.  A  homogeneous,  isotropic  sample  with  elastic  constants  appropriate 
to  basalt  was  assumed.  As  the  aspect  ratio  is  increased,  the  low  lying  modes  separate. 
For  example,  at  6  =  1.1a,  c/a  =  4,  we  expect  to  be  able  to  pick  out  fourteen  distinct 
resonances  before  mode  overlap  becomes  a  serious  problem  for  a  RUS  experiment. 

Increasing  the  aspect  ratio  further  might  allow  us  to  pick  out  even  more  dis¬ 
tinct  peaks.  However,  the  RUS  inversion  code  uses  a  fixed  order  polynomial  to  vari- 
ationally  fit  modes  [9].  As  one  side  of  a  sample  becomes  disproportionately  large,  a 
disproportionate  mimber  of  nodes  in  the  normal  modes  will  be  in  that  direction,  and 
the  inversion  code  will  lose  fitting  accuracy  in  that  direction.  We  have  chosen  to  keep 
samples  at  1/4  <  c/a  <  4  (-4  to  4  in  Fig.  3.3). 

3.5  cii  dependence 

A  rule  of  thumb  [5]  establishing  the  number  of  resonances  necessary  to  find  elastic 
constants,  is  that  five  resonance  frequencies  are  needed  to  accurately  determine  each 
independent  elastic  constant.  Thus,  for  an  isotropic  material  described  by  two  elastic 
constants,  we  need  to  experimentally  determine  the  first  ten  resonance  frequencies. 


Resonsnce  Frequency  (MHz) 
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Figure  3.3:  Calculated  resonance  frequencies  as  a  function  of  aspect  ratio.  The  sample  is 
at  a  constant  volume  of  4.8  cc;  h  =  l.ln,  and  the  aspect  ratio  is  c/a  on  the  right  side,  -a/c 
on  the  left  side.  The  elastic  constants  are  cn  =  86.6  GPa  and  C44  =  31.9  GPa. 
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Certainly  the  confidence  with  which  the  elastic  constants  can  be  determined  is  influ¬ 
enced  by  the  involvement  of  each  elastic  constant  in  the  first  ten  modes  of  the  sample. 
It  is  important  to  mention  that  the  first  ten  resonance  peaks  lie  within  a  100  kHz 
range.  Due  to  this  relatively  small  frequency  range,  the  effects  of  dispersion  have 
been  ignored. 

Since  cu  ~  2C44,  we  expect  low  frequency  modes  to  be  more  highly  dependent 
on  C44  than  cn  (in  analogy  to  the  frequencies  of  the  modes  of  a  soft  spring  versus  a  stiff 
spring  network).  Indeed,  for  a  cubic  shaped  sample  of  basalt,  the  first  eight  resonance 
frequencies  have  a  total  dependence  on  on  of  less  than  15%.  That  is,  most  low  lying 
modes  are  shear  modes,  involving  very  little  compression.  However,  the  geometry  of 
the  sample  influences  Cn  dependence  in  the  resonance  modes.  Plate-like  and  rod¬ 
like  samples  will  have  low-lying  bending  or  flexural  modes  that  are  compressional  in 
nature. 

In  Fig.  3.4,  the  cn  dependence  of  the  first  ten  resonant  modes  of  a  paral¬ 
lelepiped  sample  are  shown  as  a  function  of  the  ratio  of  the  longest  side  to  the  shortest 
side  of  the  sample  c/a.  Positive  aspect  ratios  denote  rod-like  samples  (the  aspect  ratio 
is  c/a  on  the  right  side);  negative  aspect  ratios  denote  plate-like  samples  (the  aspect 
ratio  is  —a/c  on  the  left  side).  These  dependencies  are  found  by  taking  the  deriva¬ 
tive  of  the  frequency  with  respect  to  each  of  the  moduli.  A  homogeneous,  isotropic 
sample  with  elastic  constants  appropriate  to  basalt  was  assumed.  As  the  aspect  ratio 
is  increased,  the  Cn  dependence  increases.  For  example,  for  a  sample  with  c/a  =  4, 
seven  of  the  first  ten  modes  have  a  cn  dependence  over  20%,  as  opposed  to  only  two 


Fractional  Dependence  on  c 
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Aspect  Ratio 

Figure  3.4:  cu  dependence  as  a  function  of  aspect  ratio.  The  sample  is  at  a  constant 
volume  of  4.8  cc;  b  =  l.ln,  and  the  aspect  ratio  is  c/a  on  the  right  side,  -a/c  on  the  left 
side.  The  elastic  constants  are  cn  =  86.6  Gpa  and  C44  =  31.9  GPa. 
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modes  for  c/a  =  \. 

3.6  Anisotropy 

If  isotropic  symmetry  is  broken  in  a  single  direction,  the  sample  has  hexagonal  sym¬ 
metry  and  is  called  transversely  isotropic.  Many  consolidated  materials,  such  as 
sedimentary  rock  and  laminar  systems,  are  transversely  isotropic.  Hexagonal  symme¬ 
try  is  described  by  five  independent  elastic  constants  (i.e.,  cn,  C33,  C13,  C44,  cge).  Thus, 
in  order  to  determine  hexagonal  symmetry  using  RUS,  we  might  expect  to  need  25 
resonance  frequencies.  This  is  a  prohibitively  large  number  for  low  Q  samples.  Do  we 
actually  need  this  many?  Suppose  we  have  10  measured  resonances.  Can  we  detect 
anisotropy  with  these  modes?  The  following  is  a  test  of  the  sensitivity  of  RUS  to 
anisotropy. 

Consider  the  elastic  system  described  by  the  elastic  tensor 

( ci\  C12  Ci3  0  0  0 

C\2  Cn  Ci3  0  0  0 

Cl3  <^13  C33  0  0  0 

M  =  ,  (3.3) 

0  0  0  C44  0  0 

0  0  0  0  C44  0 

V  0  0  0  0  0  ceey 
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where 


Cn  —  (!+£')  ^^0' 


ci2^0.4(l-2e)C„, 


(3.4) 


Ci3  =  0.4  (1  +  e)  Co, 

<"3.3  =  (1  ■“  26')  (7o, 

C44  =  0.3(1  -  1.5f)Co, 

C66  =  0.3  (1  +  3^)  Co, 

0  <e  <  0.25,  Co  =  10^^  dynes/cm^,  and  po  =  1  gm/cm''’.  As  £  varies  from  0  to  0.25 
the  elastic  tensor  varies  from  isotropic 


to  hexagonal 


A/(0.25)  =  Co 
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These  numbers  have  been  chosen  so  that  the  average  properties  of  the  tensor  are  the 


same  while  the  system  is  shifting  from  isotropic  to  hexagonal  with  the  change  in  s. 
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By  average,  we  mean  that  the  first  three  diagonal  components  (cn,  C22,  C33)  average  to 
Co;  the  three  ofF-diagonal  elements  (012,013,023)  average  to  0.4Co;  and  the  last  three 
diagonal  elements  (044,055,006)  average  to  0.3Co,  regardless  of  the  value  of  e.  For 
e  =  0.25  the  elastic  matrix  is  approximately  that  of  zinc  [14], 

We  proceed  as  follows: 


1.  Choose  values  of  o/a,  and  e. 

2.  Calculate  the  lowest  10  frequencies,  .91 , . . . ,  910,  using  the  elastic  tensor  in  Eqs.  (3.3) 
and  (3.4).  91, . . . ,  9io  are  the  frequencies  of  the  model  system. 

3.  Assume  the  model  system  is  isotropic  and  the  frequencies  91,..., 910  are  the 
measured  frequencies  of  the  lowest  10  modes,  the  data. 

4.  Fit  the  model  system  data  with  an  isotropic  model,  i.e.,  with  Cn  and  C44. 


In  Fig.  3.5  the  frequency  error  is  shown 


F  = 


(3.7) 


as  a  function  of  e,  where  N  =  10,  for  -4(1/4)  <  c/a  <  4,  and  c/a  is  the  ratio  of 
the  longest  (c)  to  shortest  (a)  sides.  Aspect  ratios  less  than  one  are  represented  as 
negative  reciprocals,  e.g.,  c/o.  =  1/4  is  represented  as  c/a  =  -4.  For  a  c/a  ratio  of  4 
we  see  that  the  error  due  to  attempting  an  isotropic  fit  to  data  from  an  anisotropic 
sample  remains  less  than  1%  for  e  <  0.22.  If  we  choose  1%  error  as  the  threshold 
between  a  good  fit  and  a  bad  fit,  we  do  not  have  enough  information  about  the  elastic 
properties  of  the  system  to  recognize  that  it  is  not  isotropic  if  we  are  given  only  the 
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Figure  3.5:  Frequency  error  vs.  e  for  c/a  —4,  -2,  1,  2,  3,  4. 
first  10  modes.  In  practice  the  choice  of  a  suitable  threshold  must  be  empiiically 
determined. 

The  c/a  =  2  results  in  Fig.  3.5  seems  to  imply  that  this  aspect  ratio  provides 
the  best  detection  of  anisotropy.  Figures  3.6  and  3.7  show  a  measure  of  the  sharpness 
of  the  minimum  in  F,  through  the  use  of  error  elipses,  for  two  of  the  solutions  in 
Fig.  3.5  (c/a=4,  c/a=2).  In  these  figures  the  (cu/cff ,C44/4?)  P^ir  for  which  the 
frequency  error  is  a  minimum,  is  plotted  (black  dot)  for  0  <  e  <  0.1  where  =  Co 
and  cf}  =  0.3Co.  Anisotropy  is  increasing  from  the  top  right  to  the  lower  left.  The 
black  dots  are  the  result  of  an  isotropic  fit.  The  red  ellipses  represent  a  measure  of 
the  uncertainty  in  the  isotropic  fits  (the  sharpness  of  the  corresponding  minima).  It 
follows  that  for  c/a  =  4,  the  isotropic  fit  is  very  accurately  determined  when  cn/cj? 
approaches  1,  while  the  uncertainty  grows  as  we  decrease  that  ratio  (increase  e).  From 
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Figure  3.6:  A  measure  of  the  quality  of  fit  (sharpness  of  the  minimum)  for  c/a=4  and  an 
increasing  variation  from  isotropy. 


Fig.  3.6  we  see  that  for  aspect  ratio  4  the  miniinuiii  in  the  frequency  error  is  sharply 
defined.  Contrast  this  to  Fig.  3.7  for  aspect  ratio  2.  For  this  case  the  frequency  error 
is  not  at  all  sharply  defined.  With  aspect  ratio  2  the  frequency  error  is  neither  small 
(note  it  is  less  than  1%  for  e  <  0.06)  nor  sharply  determined. 

At  e  =  0.09  (Fig.  3.8)  there  are  6  points  and  there  associated  error  ellipses 
labeled  5,  10,  15,  20,  25,  and  30.  These  points  represent  the  frequency  error  for 
c/a  =  A  found  for  isotropic  fits  of  iV  =  5, . . .  ,30  anisotropic  frequencies.  Here  we 
are  investigating  the  ability  to  see  anisotropy  from  an  isotropic  model,  as  the  number 
of  frequencies  increases.  The  differences  in  the  sharpness  of  the  minima  (size  of  the 
error  ellipses)  give  a  very  good  indication  of  how  the  ability  to  discriminate  increases 
with  an  increasing  number  of  frequencies.  We  argue  that  for  N  >  15  one  begins  to 
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Figure  3.7:  A  measure  of  the  quality  of  fit  (sharpness  of  the  minimum)  for  c/a=2  and  an 
increasing  variation  (e)  from  isotropy. 


Figure  3.8:  A  measure  of  the  quality  of  fit  (sharpness  of  the  minimum)  for  c/a=4,  e=0.09 
and  a  varying  number  of  measured  resonance  frequencies  (/  =  5,  10,  15,  20,  25,  30  as 
lableled). 
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see  errors  great,  enough  to  indicate  that  the  sample  is  not  isotropic. 

Thus  RUS  is  capable  of  distinguishing  a  deviation  from  isotropy,  given  a  suf¬ 
ficient  number  of  resonanc:e  frequencies.  A  sample  with  cj a  =  2,  displays  large  errors 
at  small  anisotropic  deviations,  as  well  as  very  shallow  minima  in  the  frequency  fits, 
thus  giving  us  an  idea  of  the  goodness  of  fit,  or  rather  that  they  don’t  fit.  For  higher 
Q  materials,  the  best  method  of  detecting  anisotropy  may  be  to  simply  measure  more 
peaks,  however,  when  the  number  of  resonance  frequencies  is  limited  to  10,  the  sample 
size  should  have  c/a  =  2  (assuming  b  =  l.ln)  to  have  the  greatest  sensitivity  to  the 


anisotropy. 


37 


Chapter  4 


Experimental  Results  And 
Conclusions 


The  previous  chapters  have  provided  a  foundation  with  which  we  can  now 
investigate  real  samples.  Using  the  analytical  method  described  in  Chap.  2,  and 
the  techniques  for  optimizing  our  results  discussed  in  Chap.  3,  we  have  performed 
RUS  experiments  on  a  variety  of  rock  types,  and  sizes.  A  case  study  of  a  sample  of 
black  gabbro  is  detailed  in  Appendix  B,  documenting  the  process  of  applying  RUS  to 
inhomogeneous  samples. 

4.1  Results 

RUS  experiments  and  analysis  were  performed  on  samples  of  berkeley  blue 
granite,  pink  quartzite,  and  black  gabbro  in  sample  set  1  (Table  4.1).  Figure  4.1 
provides  a  pictoral  representaion  of  sample  set  1.  Sample  set  2  consists  of  an  assort- 
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able  4.1:  Sample  Set  1 


Rock  Type 

ID 

a 

(cm) 

b/a 

c/a 

^/a 

Q 

(GPa) 

C44 

(GPa) 

%error 

Black  Gabbro 

BG-1 

2.8 

1.1 

1.4 

0.23 

350 

96 

36 

0.557 

Black  Gabbro 

BG-2 

2.3 

1.1 

2.6 

0.28 

350 

97 

36 

0.703 

Black  Gabbro 

BG-3 

1.0 

1.1 

4.0 

0.32 

350 

107 

35 

0.377 

Black  Gabbro 

BG-4 

1.3 

1.1 

4.0 

0.24 

350 

106 

35 

0.377 

Black  Gabbro 

BG-5 

1.6 

1.1 

4.0 

0.31 

350 

119 

35 

0.489 

Black  Gabbro 

BG-6 

2.0 

1.1 

4.0 

0.21 

350 

107 

36 

0.279 

Pink  Quartzite 

PQ-1 

2.0 

1.1 

4.0 

1.68 

250 

69.6 

35.1 

1.271 

Berkeley  Blue  Granite  BB-1 

2.0 

1.1 

4.0 

0.20 

230 

30.7 

12.3 

13.28 

ment  of  basalt,  sierra  white  granite,  and  black  gabbro.  Isotropy  was  assumed  for  all 
samples.  The  results  are  tabulated  in  order  of  increasing  aspect  ratio  c/a.  To  crudely 
determine  the  size  of  inhomogeneity  ^  in  each  sample,  the  largest  surface  structiire 
(blob)  on  each  sample  was  measured.  In  sample  set  1  cases,  the  first  ten  resonance 
frequencies  were  used  for  the  fit.  The  average  Q  is  determined  from  Q  =  u/Au;  aver¬ 
aged  over  the  first  ten  resonances,  where  cv  is  the  resonance  frequency,  and  Aa;  is  the 
full-width  at  half-maximum  of  the  resonance  intensity.  Note  that  the  RMS  percent 
error  is  strongly  dependent  on  c/a,  a/^,  and  Q.  High  values  of  these  parameters  lead 
to  better  fits  between  measured  and  calculated  resonances. 

Notice  that  for  all  black  gabbro  samples  in  sample  set  1,  the  shear  modulus 
(C44)  is  consistently  35-36  GPa.  This  agreement  between  samples  indicates  that  the 
measured  resonances  have  an  adequate  dependence  on  C44.  However,  there  is  a  marked 
difference  between  the  compressive  modulus  (cn)  for  samples  having  an  aspect  ratio 
c/o  =  4  and  those  for  which  c/a  <  4.  The  error  for  the  samples  with  c/a  <  4  is  also 
greater  than  for  the  c/a  =  4  samples.  This  indicates  that  the  lack  of  cu  dependence 
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Varying  Volume  (size) 


Figure  4.1;  Pictoral  roprosentat.ion  of  sample  set.  1 


Table  4.2:  Sample  Set  2 


Rock  Type 

ID 

a 

(cm) 

h/a 

c/a 

^/a 

Q 

<^11 

(GPa) 

C44 

(GPa) 

Sierra  White  Granite 

SW-1 

3.7 

1.1 

1.5 

0.18 

145 

37.6 

18.7 

Sierra  White  Granite 

SW-2 

2.8 

1.1 

1.6 

0.24 

140 

43.2 

18.6 

Sierra  White  Granite 

SW-3 

1.0 

1.25 

3.3 

0.46 

150 

20.8 

21.3 

Sierra  White  Granite 

SW-4 

1.7 

1.2 

4.0 

0.17 

140 

37.5 

19.1 

Basalt 

B-1 

2.1 

1.1 

1.5 

0.29 

275 

88.1 

31.3 

Basalt 

B-2 

2.7 

1.0 

1.7 

0.14 

240 

88.4 

31.7 

Basalt 

B-3 

2.8 

1.1 

2.0 

0.04 

255 

87.3 

31.6 

Basalt 

B-4 

1.7 

1.1 

4.0 

0.06 

335 

87.2 

31.5 

1.238 

2.023 

17.29 

0.630 
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in  the  low-lying  modes  for  the  more  compact  samples  is  effecting  the  final  result,  as  we 
expected  from  our  modeling  in  Sec.  3.5.  Sample  SW-2  (Table  4.2)  also  demonstrates 
this  phenomena  as  is  shown  from  the  larger  fit  error  and  the  slightly  higher  value  for 
Cii  compared  to  the  cn  value  of  samples  SW-1  and  SW-4. 

Sample  set  1  also  illustrates  the  effects  of  decreasing  Q.  The  fits  for  samples 
BB-1  and  PQ-1  (low  Q)  have  a  substantially  larger  %  error  than  the  BG  sample 
fits  (high  Q).  It  is  interesting  to  note,  however,  that  while  the  sample  BB-1  fit 
has  a  13%  error,  the  commonly  accepted  values  for  cn  and  cu  are  30  GPa  and  13 
GPa,  respectively.  Thus  even  though  the  fitting  error  is  high,  the  moduli  determined 
from  the  fit  are  quite  accurate.  It  is  important  to  remember  that  results  from  the 
RUS  inversion  are  the  best  possible  answers  for  the  given  data.  The  %  error  is  an 
uncertainty  level,  not  an  a  measure  of  the  difference  between  real  and  determined 
values  of  the  moduli. 

In  sample  set  2,  we  have  a  random  distribution  of  aspect  ratios.  Here  we 
can  see  similar  results  as  for  sample  set  1.  As  for  the  black  gabbro  samples  above, 
both  Sierra  white  granite  and  basalt  give  results  with  lower  percent  error  for  samples 
having  an  aspect  ratio  c/a  =  4.  Also  we  again  see  close  agreement  in  the  shear 
moduhis  between  samples  of  the  same  material. 

Sample  SW-3  is  an  interesting  sample  in  that  it  is  the  only  SW  sample  that 
exhibits  a  high  percent  error,  and  whose  moduli  deviate  from  the  other  SW  samples. 
This  is  not  surprising  as  we  look  at  the  size  of  the  inhomogeneity  in  the  sample.  In 
Sec.  3.2  we  stated  that  the  relative  size  of  the  inhomogeneities  (^/o)  cannot  exceed 
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1/5  (0.2).  Sample  SW-3  has  a  ^/a  value  of  0.46.  more  than  twice  that  of  what  should 
be  allowed.  In  this  sample  inhomogeneity  prevents  us  from  determining  the  average 
elastic  moduli.  The  value  for  i/a  given  in  Sec.  3.2  is  actually  a  conservative  value,  as 
the  model  used  was  one-dimensional.  In  a  three  dimensional  elastic  solid,  the  first  ten 
resonances  will  not  all  be  along  the  same  dimension,  and  thus  a  larger  inhomogeneity 
could  be  tolerated.  This  would  explain  why  the  BG  samples  all  exhibit  a  low  %  error 
yet  have  i/a  values  slightly  above  0.2. 

4.2  Conclusions 

The  modeling  in  Chap.  3  and  the  results  shown  in  Tables  4.1  and  4.2  indi¬ 
cate  that  RUS  is  a  viable  technique  for  characterizing  the  average  elastic  behavior 
of  inhomogeneous  materials.  Although  larger  RMS  errors  can  be  expected  for  inho¬ 
mogeneous  materials  than  those  acceptable  in  the  single  crystal  business  (<  0.5  % 
[5]),  our  results  are  generally  close  for  different  samples  of  the  same  material,  and 
consistent  with  accepted  values  [15]. 

For  RUS  on  rock  samples  to  be  a  viable  technique,  it  must  be  developed  further, 
specifically  to  the  needs  of  the  field  of  geophysics.  For  example,  we  need  to  develop 
RUS  for  cylindrical  sample  geometry,  so  that  it  is  capable  of  accepting  a  standard  core 
sample.  The  symmetry  of  cylinders  poses  new  problems  that  will  be  considered  in  the 
future.  Also  we  plan  to  take  a  more  in-depth  look  at  the  possibility  of  detecting  and 
analyzing  anisotropic  systems.  Finally  we  plan  to  develop  an  environmental  control 
chamber  inside  which  the  RUS  experiment  can  be  conducted.  The  chamber  will  allow 
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us  to  introduce  two  new  variables  to  the  experiment,  temperature  and  water  vapor 
pressure.  The  chamber  design  is  currenth'^  that  of  a  sealed  vacuum  oven  in  which  we 
can  easily  achieve  250°C  (423  K),  and  enclose  standard  saturation  salts. 

With  these  new  additions  to  the  work  that  has  been  completed  here,  we  expect 
to  develop  a  new  use  for  RUS  in  geophysics.  We  will  be  able  to  explore  the  elastic 
responses  of  materials  in  the  earth’s  crust  tinder  various  temperature  and  saturation 
conditions. 
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Appendix  A 


Perturbation  Theory 


The  elastic  energy  of  a  solid  body,  in  steady  state  at  frequency  u,  is  described 
by  the  Lagrangian 


L=  fdx  ^(x)  .  (A.l) 

where  u  is  the  displacement  vector  at  postition  x,  C4jifi  is  the  elastic  tensor,  p  the  mass 
density,  repeated  indices  are  summed  over  the  cartesian  coordinates,  and  0  describes 
the  extent  or  figure  of  the  sample, 


0(x) 


I  1  X  inside  the  sample 


(A.2) 


0  X  outside  the  sample 

Equation  (A.l)  is  quite  general,  allowing  for  (1)  an  arbitrary  sample  figure  0(x)  (2) 
a  nonuniform  density  p(x),  and  (3)  a  non  uniform  elastic  tensor  Cijki{x). 

The  equation  of  motion  for  the  normal  modes  is  found  by  varying  L  with 
respect  to  Uj.  If  the  traction  on  the  surfaces  defined  by  </>(x)  vanishes,  w,:  satisfies  a 
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wave  equation  in  the  form 


p(xVS„  +  (c«h(x)||)  =  0 


(A.3) 


for  X  inside  of  the  sample.  The  condition  that  the  traction  on  the  surfaces  vanish  is 
enforced  by  holding  the  sample  so  that  it  is  effectively  free  standing. 

The  equation  of  motion  for  ?/,  ,  Eq.  (A.3),  can  be  cast  in  the  form  of  a  variational 
problem  [16].  That  is,  the  quantity  where 

.1  dx  (f)(x)  Cijki(x)§^^ 


~  f  dx  (l){x)  p{x)ii^ 

must  be  stationary  subject  to  arbitrary  variations  of  u,  consistent  with  traction  free 
boundaries.  Using  this  form  for  the  normal  mode  frequencies  it  is  possible  to  make  a 
systematic  study  of  the  consequences  of  change  in  0(x),  p(x),  and  Cjjki{x).  Assume 
the  ideal  sample  is  a  rectangular  parallelepiped  specified  by  (j)o,  has  uniform  density 
po,  and  has  uniform  elastic  constants  c?  *.,.  Then  variations  in  these  quantities  are 
given  by  6(f){x)  =  <l>{x)  —  (l>o,  ^p{x)  =  p{x)  —  po  and  6cijki{x)  =  Cijki{x)  —  c^jki-  To  first 
order  in  6(j)^  6p^  and  8c  we  have 


(A.4) 


where 


2  ^0  6N^  6Dp  8D4,' 

~  Do[  ^  No  No  Do  Do  _ 

(A.5) 

.r  0  A  duiduk 

(A.6) 

Do  =  Po  f  dx  (fin  tij, 

(A.7) 

f  /  /  \  0  \  ^*6:  duk 

8Nc  —  1  dx  <f)o  \Cijki{x)  ’ 

(A.8) 

8N4,  -  c%ki  1  dx  {(f>{x)  <j)o)  , 

(A.9) 
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8 Dp  =  1  dx  00  (p(x)  -  Po)n  - . 

(A.IO) 

8D^  =  po  f  dx  {0(x)  -  0o)  vf. 

(A.ll) 

Using  the  variational  technique  of  Visscher  et  al.  we  can  find  the  eigenvalues 
uIq  =  Nq/Do  and  eigenfunctions  —  i/v(x)  associated  with  the  ideal  sample.  Thus 
the  lowest  order  contribution  to  the  frequency  shift  due  to  a  perturbation  is 


6N,  8Dp  6D4, 


(A.12) 


No  No  Do  Do 
where  u.,;  =  Vv(x).  For  the  example  of  an  inhomogeneous  mass  density  we  would  have 

(A.13) 


-  ^1,0  ^  ./  dx(^o  (p(x)  -po)|V;„(x)p 

PoJ  dx</»0 


(jlJ, 


i/,0 


where 


rO  r  dv  d»n  ^Vv(x)  gVv(x) 
2  _  ^ijkl  ■'  yO  a.r.j  dxi 

~  po.f  rfx  00  |V’t.(x)P 


(A.14) 


Equation  (A.13)  is  used  in  Sec.  3.3  for  the  case  of  a  mass  defect  to  illustrate  the 
consequences  of  a  mass  inhomogeneity  on  resonant  mode  frequencies. 
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Appendix  B 


Case  Study 


In  this  appendix  we  will  document  a  complete  experiment  from  the  forward  calculation 
to  final  result.  We  chose  to  work  with  black  gabbro  for  its  relatively  high  homogeneity 
and  Q.  Figures  B.l  and  B.2  demonstrate  why  high  Q's  are  important.  In  Fig.  B.l, 
a  sample  spectrum  for  stainless  steel  {Q  >  1000)  is  shown;  in  Fig.  B.2  a  sample 
spectrum  of  Berea  sandstone  {Q  w  50)  is  shown.  The  spectrum  in  Fig.  B.l  shows 
the  resonance  frequencies  very  clearly.  The  first  15  resonances  are  very  easy  to  pick. 
In  contrast,  the  spectrum  given  in  Fig.  B.2  contains  what  appears  to  be  a  lot  of 
noise.  Consequently  the  first  15  resonances  are  difficult  to  discern.  Most  rocks  have 
resonance  spectra  that  look  more  like  Berea  sandstone  than  stainless  steel.  In  this 
study  we  have  limited  ourselves  to  high  Q  rocks  {Q  >  150). 

The  particular  sample  of  Black  Gabbro  analyzed  in  this  appendix  is  the  fourth 
black  gabbro  sample  characterized  in  Chap.  4,  Table  4.1.  This  sample,  designated  BG- 
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4,  htvs  dimensions  1.29  x  1.40  x  5.01  c,nr\  fuspect  ratios  of  e/a  =  4  and  h/a  =  1.1,  mass 
26.18  g,  and  largest  observable  inhomogeneity  ^  =  0.47  cm.  The  dimensions,  mass, 
and  an  initial  gness  for  the  elastic  constants  (cn  =  86.6  GPa,  C44  =  31.9  GPa)  are  the 
input  for  a  numerical  forward  calculation  of  resonance  frequencies.  This  input  file  is 
shown  in  Fig.  B.3.  The  output  from  the  forward  calculation  (also  shown  in  Fig.  B.3) 
displays  the  first  15  resonance  frecpiencies  in  MHz,  and  their  fractional  dependence 
upon  the  two  independent  moduli  (df/d(moduli)  columns,  cu,  C44  respectively).  Note 
that  the  first  ten  resonances  are  well  separated  and  have  a  strong  dependence  on  the 
compressional  modulus  Cn  This  is  largely  due  to  our  choice  of  aspect  latio  c/a  —  4, 
as  suggested  by  our  numerical  studies. 

This  output  file  provides  an  initia,!  range  of  frequencies  to  guide  the  experiment. 
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Based  on  the  results  in  Fig.  B.3  the  first  frequency  scan  (Fig.  B.4)  was  performed  over 
the  range  20-100  kHz.  This  range  was  expected  to  contain  the  first  ten  resonances  and 
still  allow  for  a  shift  in  those  resonances  from  the  predicted  values.  The  first  broad 
sweep  of  frequencies  reveals  what  appears  to  be  nine  of  the  ten  frequencies  that 
the  forward  calculation  led  us  to  expect.  To  accurately  determine  the  experimental 
resonance  freqiiencies,  we  narrowed  successive  scans  to  the  regions  labeled  A-G  on 
Fig.  B.4. 

Accurate  determination  of  the  lowest  frequency  is  very  important,  as  the  first 
resonant  mode  depends  heavily  on  c\\  (see  Fig.  B.3).  Since  cn  dependence  in  resonant 
modes  is  often  difficult  to  achieve,  but  is  necessary  for  accurate  determination  of  cn, 
it  is  essential  that  we  find  the  first  resonance  of  the  sample.  Missing  peaks  can  be 
accounted  for  in  the  analysis  with  null  entries,  but  influence  the  accuracy  of  our 
experiment.  In  Fig.  B.5  (region  A  in  Fig.  B.4)  we  see  that  what  appeared  to  be  a 
single  peak  is  actually  two  closely  spaced  peaks.  Adjusting  the  phase  of  the  scan 
(seen  in  the  upper  left  window  of  each  scan)  brings  the  peaks  out  more  prominently 
for  selection.  The  selected  resonances  occur  at  24.6  kHz,  26.0  kHz.  These  values 
are  shifted  from  the  expected  values  (22.8  kHz,  24.2  kHz),  indicating  that  one  or 
more  of  the  values  entered  in  the  forward  calculation  are  incorrect.  Re-measuring 
and  weighing  the  sample,  we  found  no  error  in  the  measurements.  Consequently  we 
assume  that  our  initial  guess  of  elastic  constants  is  not  acx'urate. 

Scanning  regions  A-G  in  turn  (Figs.  B.5-  B.ll),  pinpointed  each  of  the  ex¬ 
pected  resonance  frequencies  without  confusion  due  to  artifacts  or  induced  noise.  In 
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Forward  Calculation  Input  File: 


BG-4 

2  3  10  15  26.18  1.00  1 

.866  .319 
5.01  1.29  1.4 


Resulting  Output  File: 

BG-4 

free  moduli  are  cll,  c44 

using  10  order  polynomials  mass=  26.1800  gm  rho 


n 

f  ex 

fr  %err 

wt 

k 

i 

df /d (moduli ) 

1 

1.000000 

0.022814-97.72 

1.00 

7 

2 

0.35 

0.65 

2 

1.000000 

0.024163-97.58 

1.00 

6 

2 

0.34 

0 . 66 

3 

1.000000 

0.030326-96.97 

1.00 

4 

1 

0.00 

1.00 

4 

1.000000 

0.050743-94.93 

1.00 

3 

2 

0.29 

0.71 

5 

1.000000 

0.051376-94.86 

1.00 

5 

1 

0.36 

0.64 

6 

1.000000 

0.052610-94.74 

1.00 

2 

2 

0.28 

0.72 

7 

1.000000 

0.060601-93.94 

1.00 

8 

2 

0.00 

1.00 

8 

1.000000 

0.081655-91.83 

1.00 

7 

3 

0.24 

0.76 

9 

1.000000 

0.083553-91.64 

1.00 

6 

3 

0.23 

0.77 

10 

1 . 000000 

0.090749-90.93 

1.00 

4 

2 

0.00 

1.00 

11 

1.000000 

0.101677-89.83 

1.00 

1 

2 

0.27 

0.73 

12 

1.000000 

0.109945-89.01 

1.00 

2 

3 

0.13 

0.87 

13 

1.000000 

0.110721-88.93 

1.00 

3 

3 

0.16 

0.84 

14 

1.000000 

0.121978-87.80 

1.00 

8 

3 

0.01 

0.99 

15 

1.000000 

0.131907-86.81 

1.00 

4 

3 

0.06 

0.94 

Bulk  Modulus=  0.441 

cll  c22  c33  c23  cl3  cl2  c44 

0.8660  0.8660  0.8660  0.2280  0.2280  0.2280  0.3190  0 


dl  d2  d3 

5.01000  1.29000  1.40000 


2.893  gm/cc 


c55  c66 

3190  0.3190 


Figure  B.3:  Input  file  for  sample  BG-4,  forward  calculation,  and  the  resulting  output  file. 
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Figxire  B.5;  Scan  of  region  A  for  sample  BG-4.  Decreased  scanning  range  to  resolve  the 
resonance  peak. 


Figure  B.6:  Artifact,  scan  of  sample  BG-4,  region  E.  No  resonance  peaks  in  this  range. 
Area  scanned  to  illustrate  the  obvious  difference  to  ”  peak”  shaped  artifacts  and  the  actual 
peaks  illustrated  in  the  other  scans. 


Figure  B.7:  Sample  BG-4,  scan  of  region  B. 
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Figure  B.8:  Sample  BG-4,  scan  of  region  C. 
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Figure  B.9:  Sample  BG-4,  scan  of  region  D. 


the  first  scan  (Fig.  B.4)  there  appears  to  be  a  small  deviation  from  the  background 
noise  just  before  80kHz  (region  E).  A  closer  scan  of  this  area  (Fig.  B.6)  illustrates 
the  difference  between  a  system  artifact  and  resonance.  The  bumps  in  region  E  have 
a  much  smaller  relative  amplitude  than  the  resonances  in  region  A  (Fig.  B.5),  and 
also  have  a  more  amorphous  shape  than  the  well-defined  peaks  in  region  A.  In  ad¬ 
dition,  in  subsequent  scans,  region  E  changes  shape,  while  resonance  frequencies  do 
not.  This  indicates  that  there  is  a  weak  coupling  somewhere  in  the  system,  yet  it  is 
easily  distiguishable  from  true  resonances  in  the  spectrum. 
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Figure  B.ll;  Sample  BG-4,  scan  of  region  G. 


BG-4 
2  3  10 

0.866  O 
1.29  1 . 

0.024 


O 

0 

0 

O 

0 

0 

O 

O 

0 


026 

032 

054 

055 

056 

064 

086 

089 

095 


O  26.18  1  1 

.  319 
4  5.01 

621  0 . 024762 

001  O . 028456 

160  O . 032331 

135  0 . 049300 

135  0.053805 

274  0.057833 

016  O . 064169 

995  0.085516 

145  0.088649 

878  O . 094112 


Samplo  Name 

3rcl  entry  le  polynomial  solution  order 
Sth  entry  Is  sample  mass 

Initial  Guesses  for  c11  and  c44,  respectively 
Dimensions  (order  Is  Irrelevant  for  Isotropic) 


1.00 
1 . 00 
1.00 
1.00 
1.00 
1 . 00 
1.00 
1.00 
1.00 
1.00 


l8t  &  2nd  columns  are 
experimentally  measured  frequencies 

3rd  column  Is  a  weighting 
parameter  for  each  frequency 


Figure  B.12:  Input,  file  for  sample  BG-4,  containing  experimentally  measured  frequencies 
from  figs.  B.4-  B.ll. 
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After  experimentally  measuring  the  resonance  frequencies  of  our  sample,  the 
next  step  is  to  run  the  inversion  algorithm  [8].  Figure  B.12  includes  the  same  sample 
information  as  the  input  file  used  for  the  forward  calculation  (Fig.  B.3)  with  the 
addition  of  three  columns  containing  the  measured  frequencies  (columns  one  and 
two)  and  a  weighting  column  indicating  the  confidence  level  of  the  measurement 
(0  <  weighting  value  <  1).  The  final  output  is  shown  in  Fig.  B.13.  The  information 
provided  in  the  output  file  includes  the  measured  resonance  frequencies,  expected 
resonance  frequencies,  the  difference  between  them  (fex,  fr,  %err  columns),  the  bulk 
modulus  (an  engineering  constant),  individual  moduli,  and  the  accuracy  of  the  final 
solution  (rms  error). 
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BG-4 

free  moduli  are  cll,  c44 

free  dimensions  are  dl,  d2 ,  d3 

using  10  order  polynomials  mass=  26.1800  gm  rho=  2.893  gm/cc 


n 

f  ex 

fr 

%err 

wt 

k 

i 

df/d (moduli)  | 

1 

0 . 024621 

0.024533 

-0.36 

1. 

1 

2 

0.28 

0.72  i 

2 

0 . 026001 

0.025976 

-0.10 

1. 

7 

2 

0.28 

0.72 

3 

0 . 032160 

0 . 032055 

-0.33 

1. 

4 

1 

0.00 

1.00  ! 

4 

0 . 054135 

0 . 054361 

0.42 

1. 

2 

2 

0.23 

0.77 

5 

0 . 055135 

0 . 055040 

-0.17 

1. 

5 

1 

0.29 

0.71 

6 

0 . 056274 

0 . 056326 

0.09 

1. 

8 

2 

0.22 

0.78 

7 

0 . 064016 

0.064057 

0.06 

1. 

3 

2 

0.00 

1.00  i 

8 

0 . 086995 

0 . 087226 

0.27 

1. 

1 

3 

0 . 19 

0.81  ; 

9 

0 . 089145 

0 . 089194 

0.05 

1 . 

7 

3 

0.18 

0.82 

10 

0 . 095878 

0.095932 

0 .06 

1. 

4 

2 

0.00 

1.00 

Bulk  Modulus=  0.573 

cll  c22  c33  c23  cl3  cl2  c44  c55  c66 

1.0444  1.0444  1.0444  0.3370  0.3370  0.3370  0.3537  0.3537  0.3537 


dl  d2  d3 

1.29241  1.40270  4.99105 


loopttlO  rms  error=  0.2319  %,  changed  by  -0.0000063  % 
length  of  gradient  vector=  0.194379  blamb=  0.000000 


eigenvalues 
0 . 02097 
0 . 81465 
2.66530 
34.43162 
9827496.70739 


eigenvectors 

0.97  0.25-0.04-0.06  0.00 
0.05-0.43  0.00-0.90  0.00 
-0.03  0.20  0.65-0.10-0.72 
-0.04  0.02-0.75-0.01-0.67 
0.24-0.85  0.13  0.41-0.19 


chisquare  increased  2%  by  the  following  %  changes  in 
independent  parameters 

0.94  0.15  -0.05  -0.05  0.10 
0.04  -0.20  0.05  0.00  -0.06 
0.00  0.00  0.09  -0.10  0.00 


Figure  B.13:  Final  output  file  for  sample  BG-4. 
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